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Abstract 

We calculate the density-dependent spin-isospin asymmetry energy J(kf) of nuclear 
matter in the three-loop approximation of chiral perturbation theory. The interaction 
contributions to J(kf) originate from one-pion exchange, iterated one-pion exchange, and 
irreducible two-pion exchange with no, single, and double virtual A-isobar excitation. We 
find that the approximation to l7r-exchange and iterated l-zr-exchange terms (which leads 
! already to a good nuclear matter equation of state by adjusting an emerging contact- 

term) is spin-isospin stable, since J(fc/o) — 24MeV > 0. The inclusion of the chiral irNA- 
| ^ . dynamics, necessary in order to guarantee the spin-stability of nuclear matter, keeps this 

property intact. The corresponding spin-isospin asymmetry energy J(kf) stays positive 
' even for extreme values of an undetermined short-distance parameter J5 (whose possible 

range we estimate from realistic NN-potentials). The largest positive contribution to 
J(kf) (a term linear in density) comes from a two-body contact-term with its strength 
{SJ ' fitted to the empirical nuclear matter saturation point. 
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In recent years a novel approach to the nuclear matter problem has emerged. Its key element 
is a separation of long- and short-distance dynamics and an ordering scheme in powers of small 
momenta. At nuclear matter saturation density po — 0.16fm -3 the Fermi momentum kf and 
the pion mass m n are comparable scales (fc/o — 2m n ), and therefore pions must be included 
as explicit degrees of freedom in the description of the nuclear many-body dynamics. The 
contributions to the energy per particle E(kf) of isospin-symmetric (spin-saturated) nuclear 
matter as they originate from chiral pion-nucleon dynamics have been computed up to three- 
loop order in Refs.fl, 2]. Both calculations are able to reproduce the empirical saturation point 
of nuclear matter by adjusting one single parameter (either a contact-coupling go + gi — 3.23 [1] 
or a cutoff scale A ~ 0.65 GeV [2]) related to unresolved short- distance dynamics. 1 The basic 
mechanism for saturation in this approach is a repulsive contribution to the energy per particle 
E(kf) generated by Pauli-blo eking in second order (iterated) pion-exchange. As outlined in Sec. 
2.5 of Ref.[2] this mechanism becomes particularly transparent by taking the chiral limit m n = 0. 
In that case the interaction contributions to E(kf) are completely summarized by an attractive 
fc^-term and a repulsive fci-term where the parameter-free prediction for the coefficient of the 
latter is very close to the one extracted from a realistic nuclear matter equation of state. 

In a recent work [3] we have extended the chiral approach to nuclear matter by including 
systematically the effects from 27r-exchange with virtual A(1232)-isobar excitation. The physi- 
cal motivation for such an extension is threefold. First, the spin-isospin-3/2 A(1232)-resonance 

1 Fitting a cutoff scale, as done in Ref.[2], must be viewed as a short-term intermediate step before an eventual 
full effective field theory calculation. Cutoff independence of physical observables is in fact a primary goal of 
effective field theory. 
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is the most prominent feature of low-energy 7riV-scattering. Secondly, it is well known that 2tt- 
exchange between nucleons with excitation of virtual A-isobars generates the needed isoscalar 
central NN-attraction [4] which in phenomenological one-boson exchange models is often sim- 
ulated by a fictitious scalar "<t" -meson exchange. Thirdly, the delta- nucleon mass splitting 
A = 293 MeV is of the same size as the Fermi momentum fc/o — ^m^ at nuclear matter satu- 
ration density and therefore pions and A-isobars should both be treated as explicit degrees of 
freedom. A large variety of nuclear matter properties has been investigated in this extended 
framework in Ref.[3]. It has been found that the inclusion of the chiral 7riVA-dynamics is able 
to remove most of the shortcomings of previous chiral calculations of nuclear matter [2, 5, 6, 7]. 
However, there remain open questions concerning the role of yet higher orders in the small mo- 
mentum expansion and its "convergence". The relation of the fitted short- distance parameters 
[2] to those of few-nucleon systems is not clear at this moment. Also, a rigorous power counting 
that justifies the perturbative chiral expansion for nuclear matter has not yet been formulated. 
Recent work by Bogner et al. [8] based on the universal low-momentum NN-potential Vj ow -k 
may open interesting perspectives in this direction. 

Irrespective of such foundational questions it is also necessary to check various stability 
conditions for nuclear matter in the chiral framework. In a recent paper [9] we have ana- 
lyzed spin-stability. It turned that the inclusion of the chiral 7riVA-dynamics is essential in 
order to guarantee the spin-stability of isospin-symmetric nuclear matter. The truncation to 
fourth order terms in the small momentum expansion with interaction contributions only from 
In- and iterated l7r-exchange is spin-unstable [9]. This statement holds independently of the 
regularization scheme if the contact-terms (generating contributions linear in the nucleon den- 
sity) are consistent with the empirical nuclear matter bulk properties: E(kf ) ~ —16 MeV and 
A(kfo) ~ 34 MeV. Now, since a nucleon possesses four internal spin and isospin degrees of 
freedom one can prepare nuclear matter also in a spin-isospin mixed asymmetric configuration. 
The stability of nuclear matter against such correlated spin-isospin deformations is the sub- 
ject of the present paper. For recent work on generalized symmetry energy coefficients in the 
context of phenomenological Skyrme forces, also see Ref.flO]. Analogous earlier studies within 
Brueckner theory using the Reid soft-core NN-potential can be found in Ref.[ll] . 

Let us begin with defining the spin-isospin asymmetry energy J(kf) of (infinite) nuclear 
matter. Consider this many-nucleon system in a state where the equal densities of the spin- 
up protons (pt) an d the spin-down neutrons (nj) have an excess over the equal densities of 
the spin-down protons (pi) and the spin-up neutrons (nf). With the help of the spin- and 
isospin-projection operators, (1 ± cr 3 )/2 and (1 ± r 3 )/2, such a spin-isospin mixed asymmetric 
configuration is realized by the substitution: 

o(k f -\p\) - l±™ e (k + -\p\) + ±^e{k--\p\), (l) 

in the medium insertion. 2 Here, k + = kf(l + e) 1 / 3 and A;_ = kf(l — e) 1 / 3 (with e a small 
parameter) are different Fermi momenta, chosen such that the total nucleon density p = (/c 3 + 
/c 3 )/37r 2 = 2k^/3n 2 stays constant. Note that Eq.(l) describes a rather peculiar asymmetric 
configuration of nuclear matter with equal densities of protons, neutrons, spin-up states and 
spin-down states: p v = p n = P\ = P\ = k^f&ix 1 . The expansion of the energy per particle of 
spin-isospin polarized nuclear matter: 

E(k + , fc_) gr _ pol = Ejkj) + e 2 J(k f ) + 0(e 4 ) , k± = k f (l± e) 1 / 3 , (2) 

2 Medium insertion is a technical notation for the difference between the in-medium and vacuum nucleon 
propagator [2]. Effectively, it sums hole-propagation and the absence of particle-propagation below the Fermi 
surface \p\ < kf. 
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Figure 1: The two-loop one-pion exchange Fock diagram and the three- loop iterated one-pion 
exchange Hartree and Fock diagrams. The combinatoric factors of these diagrams are 1/2, 1/4 
and 1/4, in the order shown. 



defines the spin-isospin asymmetry energy J(kf). The obvious criterion for the spin-isospin 
stability of nuclear matter is then the positivity of the spin-isospin asymmetry energy: J(kf) > 
0. The energy per particle at fixed nucleon density p must take on its absolute minimum value 
in the spin- and isospin-saturated configuration. 

The first contribution to the spin-isospin asymmetry energy J{kf) comes from the kinetic 
energy y/M 2 + p 2 — M of a non-interacting relativistic Fermi gas of nucleons: 



J(k f ) 



6M 12M 3 



k f 



(3) 



with M = 939 MeV the (average) nucleon mass. The next term in this series, kj/16M 5 , is 
negligibly small at the densities of interest. 

Next, we come to interaction contributions to J(kf). The closed in-medium diagrams re- 
lated to one-pion exchange (Fock diagram) and iterated one-pion exchange (Hartree and Fock 
diagrams) are shown in Fig. 1. Differences in comparison to the calculation of the energy per 
particle E(kf) in Ref.[2] occur only with respect to the factors emerging from the spin and 
isospin traces over closed nucleon lines and the radii k± = kf(l ± e) 1 / 3 of the Fermi spheres to 
be integrated over. After some analytical calculation we find the following contribution to the 
spin-isospin asymmetry energy J(kf) from the l7r-exchange Fock diagram in Fig. 1 (including 
its relativistic l/M 2 -correction): 
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Here, we have introduced the abbreviation u = kf/m n where = 135 MeV stands for the 
(neutral) pion mass. As usual f n = 92.4 MeV denotes the weak pion decay constant and we 
choose the value qa = 1-3 of the nucleon axial-vector coupling constant in order to have a pion- 
nucleon coupling constant of g n N = QaM/ f w = 13.2. In the second and third diagram in Fig. 1 
the l7r-exchange interaction is iterated (once) with itself. These second order diagrams carry 
the large scale enhancement factor M (the nucleon mass). It stems from an energy denominator 
that is equal to a difference of small nucleon kinetic energies. With a medium insertion at each 
of two equally oriented nucleon propagators we obtain from the three-loop Hartree diagram in 
Fig. 1 the following contribution to the spin-isospin asymmetry energy: 



J(k f ) 



irg A Mml 
6(4tt/ w ) 4 



15m + -^j ln(l + 4m 2 ) - 14m - 16m 2 arctan 2m J 



(5) 
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The right Fock diagram of iterated l7r-exchange (see Fig. 1) with two medium insertions on 
non-neighboring nucleon propagators gives rise on the other hand to a contribution to the 
spin-isospin asymmetry energy of the form: 



J(k f ) = 



ng%Mm*(21u 
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This expression does not include the contribution of a linear divergence f£° dl 1 of the momentum- 
space loop integral. In dimensional regularization such a linear divergence is set to zero, whereas 
in cut-off regularization it is equal to a momentum space cut-off A. The additional term spe- 
cific for cut-off regularization will be given in Eq.(13). An in-medium diagram with three 
medium insertions represents Pauli-blocking effects in intermediate NN-states induced by the 
filled Fermi sea of nucleons. The unequal filling of the ip\,n\) and {pl,n\) Fermi seas shows 
its consequences in the spin-isospin asymmetry energy. After some extensive algebraic manip- 
ulations we end up with the following double-integral representation of the contribution to the 
spin-isospin asymmetry energy J{kf) from the Hartree diagram in Fig. 1 with three medium 
insertions: 



J(k f ) 
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where we have introduced several auxiliary functions: 

s = xy + \Ju 2 — x 2 + x 2 y 2 , ^ 
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Note that Eq.(7) stems from a nine-dimensional principal- value integral over the product of 
three Fermi spheres of varying radii k± = kf(l ± e) 1 / 3 which has been differentiated twice 
with respect to e at e = 0. Of similar structure is the contribution to J(kf) from the iterated 
l7r-exchange Fock diagram with three medium insertions. Because of the two different pion 
propagators in the Fock diagram one ends up (partially) with a triple-integral representation 
for its contribution to the spin-isospin asymmetry energy: 
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Here, we have split into factorizable and non-factorizable parts. These two pieces are distin- 
guished by whether the (remaining) nucleon propagator in the three-loop Fock diagram can be 
canceled or not by terms from the product of 7riV-interaction vertices. The factorizable terms 
can be expressed through the auxiliary function: 

G = u(l + u 2 + x 2 ) - J_[l + ( u + xf] [1 + (u - xf] In ] + ^ + X \l , (10) 

Ax 1 + (u — x) A 

and its partial derivatives for which we have introduced a (short-hand) double-index notation: 

• • d i+j G 

Gn = xV - - . , Ki + j<2. (11) 

For the presentation of the nonfactorizable terms one needs also copies of the quantities s and 
s' defined in Eq.(8) which depend (instead of y) on another directional cosine z: 

dt 



t = xz + Vu 2 — x 2 + x 2 z 2 , t' — u — . (12) 

ou 

In the chiral limit m n = the fourth order contributions in Eqs.(5-9) sum up to a negative 
fc)-term of the form: J(k f )\ mn=0 = -(^A;//47r^) 4 (M/405)(327r 2 + 741 + 1848 ln2). Finally, we 
give the expression for the linear divergence specific to cut-off regularization: 



to which only the iterated l7r-exchange Fock diagram (with two medium insertions) has con- 
tributed. In the case of the Hartree diagram the linear divergence drops out after taking the 
second derivative with respect to e. One observes that the term in Eq.(13) is just —1/3 of the 
corresponding contribution to the energy per particle E(kf) (see Eq.(15) in Ref.[2]). In this con- 
text it is interesting to note that for terms linear in density p the relation 3J(kf)n n = —E(kf)n n 
holds generally. It is a consequence of the spin-isospin structure 3 — U\ • <? 2 f\ • r 2 of a Fierz- 
antisymmetric NN-contact interaction (see e.g. Eq.(35) in Ref.[12]). 

Now we can turn to numerical results. In Fig. 2 we show the spin-isospin asymmetry energy 
J(kf) of nuclear matter as a function the nucleon density p = 2k^/3ir 2 . The solid line corre- 
sponds to a calculation up to fourth order in small momenta. It includes besides the kinetic 
energy term Eq.(3) the contributions from static l7r-exchange and iterated l7r-exchange. For 
reasons of consistency we have dropped the small relativistic l/M 2 -correction in Eq.(4) since 
it is of fifth order in the small momenta kf and m n . The cut-off scale A = 0.61 GeV has been 
adjusted 3 to the nuclear matter saturation point po = 0.173 fm -3 and E(kfo) = — 15.3MeV. 
The resulting value of the nuclear matter compressibility K = k 2 E"(kfo) = 252 MeV is con- 
sistent with a recent extrapolation from giant monopole resonances of heavy nuclei [13], which 
gave K = (260 ±10) MeV. One can read off from Fig. 2 a positive value of the spin-isospin 
asymmetry energy at saturation density: J(kf ) = J(2m n ) = 23.9 MeV. It indicates the spin- 
isospin stability of nuclear matter in this approximation. The largest positive contribution to 
Ji^m^) = 23.9 MeV comes from the term, Eq.(13), linear in density and amounts to 59.2 MeV at 
saturation density kfo = 2m 7T . Compared to that the largest negative contribution is —32.8 MeV 
and it stems from the iterated l7r-exchange Fock diagram with two medium insertions, Eq.(6). 
The remaining numerically smaller contributions cancel each other to a large extent. It must 

3 This provisional procedure introduces a model-dependence that lies outside effective field theory. 
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Figure 2: The spin-isospin asymmetry energy J(kf) of nuclear matter versus the nucleon density 
p = 2k 3 f/3ir 2 . The solid line shows the result of a calculation up to fourth order in small 
momenta including l7r-exchange and iterated l7r-exchange. The cut-off scale A = 0.61 GeV has 
been adjusted to the saturation point: po — 0.173 fm~ 3 , E(kfo) = — 15.3MeV. The positive 
values of J(kf) indicate the spin-isospin stability of nuclear matter in this approximation. 



however be stressed that at this level of approximation, with interaction terms only from lir- 
exchange and iterated l7r-exchange, nuclear matter is spin-unstable [9]. The inclusion of higher 
order terms (in particular 27r-exchange with virtual A-isobar excitation) is mandatory in order 
to achieve spin-stability of nuclear matter. 

Therefore, we turn now to contributions to J(kf) of fifth order in the small momentum 
expansion. At three-loop order these terms are generated by (irreducible) two-pion exchange 
between nucleons. The corresponding one-loop diagrams for elastic NN-scattering are shown 
in Fig. 3. Since we are counting the delta-nucleon mass splitting A = 293 MeV (together with 
kf and small momentum scale the diagrams with single and double virtual A(1232)- 

isobar excitation shown in Fig. 4 belong to the same order. By closing the two open nucleon 
lines of the one-loop diagrams in Figs. 3,4 to either two or one ring one gets (in diagrammatic 
representation) the Hartree or Fock contribution to the energy density. The Hartree contri- 
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Figure 3: One- loop diagrams of irreducible two-pion exchange between nucleons. 
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Figure 4: One-loop two-pion exchange diagrams with single and double A(1232)-isobar excita- 
tion. Diagrams for which the role of both nucleons is interchanged are not shown. 

bution to the spin-isospin asymmetry energy J(kf) vanishes identically because the relevant 
27r-exchange NN T-matrix in forward direction is spin- independent [4, 12]. The Fock contri- 
bution on the other hand is obtained by integrating the spin- and isospin-contracted T-matrix 
over the product of two Fermi spheres of radii k± = kf(l ± e) 1//3 . We separate regularization 
dependent short-range contributions to the T-matrix (originating from the ultra-violet diver- 
gences of the one-loop diagrams in Figs. 3,4) from the unique long-range terms with the help 
of a twice-subtracted dispersion relation. The occurring subtraction constants give rise to a 
contribution to the spin-isospin asymmetry energy of the form: 



k 3 k 5 



(14) 



The dimensionless parameters B 3 = —7.99 has been adjusted in Ref.[3] to the saturation 
minimum E(kf ) = — 16MeV. Again, we recognize in the first part of Eq.(14) the relation 
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constant J5 in front of the /c^/M 4 -term is (a priori) not constrained by any empirical (ground- 
state) property of nuclear matter. The long-range parts of the 27r-exchange (two-body) Fock 
diagrams can be expressed as a dispersion-integral: 
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where ImV^, ImW 7 ^, ImV^ and ImWV are the spectral functions of the isoscalar and isovector 
central and tensor NN-amplitudes, respectively. Explicit expressions of these imaginary parts 
for the contributions of the triangle diagram with single A-excitation and the box diagrams 
with single and double A-excitation can be easily constructed from the analytical formulas 
given in Sec. 3 of Ref.[4]. The pi- and /cj-dependent weighting functions in Eq.(15) take care 
that at low and moderate densities this spectral-integral is dominated by low invariant 7r7r- 
masses 2m n < \i < 1 GeV. The contributions to the spin-isospin asymmetry energy J(kf) from 
irreducible 27r-exchange (with only nucleon intermediate states, see Fig. 3) can also be cast into 
the form Eq.(15). The corresponding non- vanishing spectral functions read [12]: 
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Figure 5: Hartree and Fock three-body diagrams related to 27r-exchange with single virtual 
A-isobar excitation. They represent interactions between three nucleons in the Fermi sea. The 
combinatoric factor is 1 for each diagram. 



ImVr(z/f) 



(17) 



Next, we come to the additional 27r-exchange three-body terms which arise from Pauli 
blocking of intermediate nucleon states (i.e. from the (1 ± cr 3 T 3 )6(k± — \p\) terms in the in- 
medium nucleon propagators [2]). The corresponding closed Hartree and Fock diagrams with 
single virtual A-excitation are shown in Fig. 5. The contribution of the left three-body Hartree 
diagram to the spin-isospin asymmetry energy J(kf) has the following analytical form: 
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The delta propagator shows up in this expression merely via the (reciprocal) mass-splitting 
A = 293 MeV. Furthermore, we have already inserted in Eq.(18) the empirically well-satisfied 
relation g^NA = ^9wn/V^ for the 7riVA-coupling constant. The parameter ( = —3/4 has 
been introduced in Sec. 2 of Ref. [3] in order to reduce a too strongly repulsive p 2 -term in 
the energy particle E(kf). It controls the strength of a three-nucleon contact interaction ~ 
((g\/Af£) (NN) 3 which has the property that it contributes equally but with opposite sign to 
the energy per particle E{kf) and the spin-isospin asymmetry energy J{kf). The contribution 
of both three-body Fock diagrams in Fig. 5 to the spin-isospin asymmetry energy J(kf) can be 
represented as: 
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Figure 6: The spin-isospin asymmetry energy J(kf) of nuclear matter versus the nucleon density 
p = 2kj-/3ir 2 . In comparison to Fig. 2 the effects from 27r-exchange with single and double virtual 
A-isobar excitation are now included. The solid, dashed, and dashed-dotted curves correspond 
to the choices J5 = 0, —0.95 and —19 of the short-range parameter J 5 introduced in Eq.(14). 
The positive values of J(kf) ensure the spin-isospin stability of nuclear matter. 

The double-indices on Gs and Gt have the same meaning as explained in Eq.(ll) for the 
function G. 

In Fig. 6 we show again the spin-isospin asymmetry energy J(kf) of nuclear matter as a 
function of the nucleon density p = 2k^/3n 2 . The solid line includes all the contributions from 
chiral In- and 27r-exchange written down in Eqs. (3-9, 14-19). The (yet undetermined) short- 
range parameter J 5 has been set to zero, J5 = 0. We note as an aside that the term linear in the 
density and the cut-off A, Eq.(13), is now not counted extra since the parameter B3 = —7.99 
[3] collects all such possible terms. Numerically, these two terms linear in density are anyhow 
almost identical. One observes in Fig. 6 a positive spin-isospin asymmetry energy J(kf) which 
rises monotonically with the density p. The inclusion of the chiral 7riVA-dynamics does therefore 
not disturb the spin-isospin stability of nuclear matter. It is also interesting to look at numerical 
values of J{kf) and their decomposition. At a Fermi momentum of kf — 2m n (corresponding to 
p = 0.173 fm -3 ) the spin-isospin asymmetry energy is now J(2m n ) = 69.5 MeV (setting J5 = 0). 
The most significant changes in comparison to the previous fourth order calculation come from 
the two-body Fock and three-body Hartree contributions Eqs. (15,18) which amount together 
to 30.6 MeV + 20.7 MeV = 51.3 MeV. About one third thereof (namely 16.6 MeV) stems from 
the three-body contact interaction proportional to ( = —3/4. 

The size of the short- distance parameter J 5 in Eq.(14) is still open and large negative values 
could endanger the spin-isospin stability. In order to get an estimate of J 5 we bring into play the 
complete set of four-nucleon contact-couplings written down in Eqs. (3,4) of Ref. [14]. This set 
represents the most general short-range NN-interaction quadratic in momenta and it involves 
seven low-energy constants Ci,...,Cj. After computing the spin-isospin asymmetry energy 
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J(kf) from the corresponding contact-potential in Hartree-Fock approximation we find: 
Af 4 Af 4 

Js = IS^ 2 " 4C<l) = 144^[ 3C,(1Pl) + C(3Po) + 3C<(3Pl) + 5C<(3P2) 1 • (22) 
In the second line of Eq.(22) we have reexpressed the relevant linear combination of Ci j2 through 
the so-called spectroscopic low-energy constants which characterize the short-range part of the 
NN-potential in the spin-singlet and spin-triplet S- and P-wave states. In that representation 
we obtain from the entries of table IV in Ref. [14] for the three NN-potentials 4 CD-Bonn, 
Nijm-II, and AV-18 the numbers: J 5 = —1.34, —0.57, and —0.94. The dashed line in Fig. 6 
shows the spin-isospin asymmetry energy J(kf) which results from taking their average value 
J5 = —0.95. The corresponding reduction of the spin-isospin asymmetry energy is negligible. 
The dashed-dotted curve in Fig. 6 corresponds to the extreme choice J5 = —19. One can see 
that even with such a large negative Js-value the spin-isospin stability of nuclear matter remains 
still preserved. We can therefore conclude that spin-isospin stability is a robust property of 
the chiral approach to nuclear matter (at least in the three- loop approximation). This is a 
important finding. 

In summary we have investigated in this work the spin-isospin stability of nuclear mat- 
ter in the framework of chiral perturbation theory. For that purpose we have calculated the 
density-dependent spin-isospin asymmetry energy J(kf) of nuclear matter to three-loop order. 
The interaction contributions to J{kf) originate from l7r-exchange, iterated l7r-exchange, and 
(irreducible) 27r-exchange with no, single, and double virtual A-isobar excitation. We have 
found that the approximation to l7r- and iterated l7r-exchange terms is spin-isospin stable, 
since J(kfo) > 0. The inclusion of the chiral 7riVA-dynamics (necessary to ensure the spin- 
stability [9] of nuclear matter) keeps this property intact. The largest positive contribution 
to J(kf) comes from a two-body contact interaction with its strength fitted to the empirical 
nuclear matter saturation point. 
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